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Based on a microscopic system reservoir model,where the associated bath is not in thermal equi- 
librium, we simulate the nonstationary Langevin dynamics and obtain the generalized nonstationary 
fluctuation dissipation relation (FDR) which asymptotically reduces to the traditional form. Our 
00 j Langevin dynamics incorporates non-Markovian process also, the origin of which lies in the decaying 

term of the nonstationary FDR. We then follow the stochastic dynamics of the Langevin particle 
f^*) ■ based on the Fokker-Planck-Smoluchowski description, in ratchet potential to obtain the steady and 

£SJ ' time dependent current in an analytic form. We also examine the influence of initial excitation 

and subsequent relaxation of bath modes on the transport of the Langevin particle to show that 
the nonequilibrium nature of the bath leads to both strong non-exponential dynamics as well as 
nonstationary current. 
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INTRODUCTION 



o 

The traditional theory of noise induced transport deals with a Langevin equation describing the motion of a 
model Brownian particle in an external periodic potential, spatially symmetric or asymmetric The nature of 
asymmetry of the external force field, in which the Brownian particle is moving, is crucial in generating biased directed 
motion. While moving in a symmetric potential, the Brownian particle is unable to generate motion in a preferred 
direction due to the detailed balance principle which can be broken easily by applying an external time dependent 
perturbation, either deterministic or random. The correlation time of the external perturbation needs to be greater 
than the correlation time of the fluctuations which the system experiences from its immediate surroundings, the 
i heat bath. A general approach in this direction involves the application of a time periodic deterministic field or the 

application of a colored noise to the system of interesti^^. Adopting a different approach, one can create directed 
C*~) ■ motion by putting the Brownian particle in a biased asymmetric periodic potential from the very beginning. The 
spatial bias in the potential is able to overcome the detailed balance principle and hence can generate motion in a 
preferred directional 3 .!^ 

The theory of directed motion has gained wide interdisciplinary attention to model the phenomena of noise in- 
duced transport, where the interplay of fluctuations and non-linearity of the system plays an important ro\e^2£&£. 
Exploitation of the nonequilibrium fluctuations present in the medium helps to generate directed motion of the Brow- 
nian particle. Presence of spatial anisotropy in the potential together with nonequilibrium perturbation enables one 
to extract useful work from random fluctuations without violating the second law of thermodynamics!^ This leads 
to its wide applicability in explaining the mechanism of molecular motors^ 1 ^ tunneling in a Josephson junction^ 
rotation of dipoles in a constant field? 3 - phase locked loop,- directed transport in photovoltaic and photoreflective 
^ ' materials^ and the efficiency of tiny molecular machine in a highly stochastic environment Motor proteins like 
kinesins,dyenins and myosins are versatile biomolecular shuttle cargo encapsulated in vesicles and are present in the 
different parts of the cell. In living cells, transport occurs via the cytoskeletal filaments and motor proteins^ Motor 
proteins are also important ingredients of the mechanism of muscle contraction and cell division^ The search for 
physical principles that enable such tiny molecular machines to function efficiently in a highly Brownian regime and 
construction of artificial molecular rotors which produce controlled directional motion mimicking molecular motor 
proteins^ are the subject of ongoing interest. 

During the last two decades, several theoretical models have been proposed using the idea of a Brownian particle 
moving in a ratchet potentiali^i 3 ^ to explain the transport mechanism under various nonequilibrium situations. The 
ratchet model and its many variants like rocking ratchet, 1 diffusion ratchet^ correlation ratchet^ flashing ratchet^ 
etc., have found wide attention in recent days£ To get a unidirectional current, either spatially asymmetric periodic 
potentials or time asymmetric external forces are necessary in these models. In explaining the above mentioned 
directional transport phenomena, most of the theoretical approaches adopt phenomenological models. The first self 
consistent microscopic attempt was made by Millonas^ in the context of construction of a Maxwell's demon like 
information engine that extracts work from a heat bath. In this microscopic construction, the Hamiltonian for the 
whole system includes a subsystem, a thermal bath and a nonequilibrium bath that represents an information source 
or sinkji^ 
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In this article, we consider a simple variant of the system reservoir hamiltonianis to model the directional transport 
processes where the associated bath is in a nonequilibrium state. The model incorporates some of the features of 
Langevin dynamics with a fluctuating barrier^ and the kinetics due to space dependent friction along with the 
presence of local hot spots j 16 ' 17 i 18 Since the theories of transport processes traditionally deal with stationary bath, 
the nonstationary transport processes have remained largely overlooked so far. We specifically address this issue and 
examine the influence of initial excitation and subsequent relaxation of bath mode o 14 ' 19 ' 20 ' 21 ' 22 ' 23 on the transport of 
system particle. We show that relaxation of the nonequilibrium bath modes may result in strong non-exponential 
kinetics and nonstationary current. The physical situation that has been addressed is that at t = 0_, the time just 
before the system and the bath are subjected to an external excitation, the system is appropriately thermalized. 
At t = 0, the excitation is switched on and the bath is thrown into a nonstationary state which behaves as a 
nonequilibrium reservoir. We follow the stochastic dynamics of the system mode after t > 0. The separation of the 
time scales of the fluctuations of the nonequilibrium bath and the thermal bath to which it relaxes, is such that the 
former effectively remains stationary on the fast correlation of the thermal noised 

The organization of the paper is as follows: We discuss in Sec. II a microscopic model necessary to compute the 
transient transport process where the system in question is not initially thermalized and the associated bath is 
thrown into a nonequilibrium and nonstationary situation by sudden initial excitation of some of the bath modes. 
Appropriate elimination of the reservoir degrees of freedom leads to a non-Markovian Langevin equation, stochasticity 
being contributed by both additive thermal noise and the multiplicative noise due to relaxing nonequilibrium modes. 
In Sec. Ill, following the prescription of Refl24l. the Fokker-Planck description is provided in position space which is 
valid for state dependent dissipation. We then derive the time dependent solution of the associated Smoluchowski 
equation for probability density function. As an application of our development, in Sec. IV, we consider the motion 
of a Langevin particle in a periodic ratchet potential and obtain the stationary and time dependent average velocity 
of the Langevin particle and show that for symmetric periodic potential, the direction of average velocity depends on 
the initial excitation of intermediate bath modes. Summarizing remarks are presented in Sec.V. 

II. THE BACKGROUND AND THE MODEL 

To make the paper self contained, we first discuss the essential features of the traditional theory of system reservoir 
dynamics in this section and then describe the model we adopt in the present work. This shows how our model 
deviates from the usual system reservoir theory and brings up the new features of our model. 

A. The traditional system reservoir model 

In the traditional system reservoir model ) 25 ' 26 ! 27 the reservoir is assumed to be in equilibrium at t = in the presence 
of the system, and the appropriate distribution of the initial state of the heat bath is governed by the Hamiltonian 

Hb + Hsb = /"] 

which includes the static interaction part, Hsb, between the system and the reservoir. The total Hamiltonian of the 
system plus bath is then usually written as 

9 

H = V — + V{x)+H B + H SB . (2) 

In Eqs.lJTJSJl, the system (mass weighted) is described by the coordinate x and the conjugate momentum p, and the 
heat bath, composed of a set of linear harmonic oscillators, by the coordinate q v and the conjugate momenta p„, 
v = 1, 2 • • • N. m v is the mass of the v-th oscillator and w„, the corresponding frequency. The system bath interaction 
is generally taken to be linear in both the system and the bath coordinates through the coupling constant g v . V(x) 
represents the external force field in which the Brownian particle is executing random motion. The bath is assumed 
to be in thermal equilibrium at temperature T and the initial distribution is considered to be a canonical on o 25 ' 26 ' 27 

W[q(0),p(0)] = |exp (- HB *T SB ) > ( 3 ) 

where, Z is the normalization constant and k B is the Boltzmann constant. To derive the dynamical equations for the 
system in terms of x and p, one usually eliminates the bath degrees of freedom from the equations of motion of the 



'l_ m v ul f _ g v x \ 
n v 2 \ m v Jl ) 
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system variable 2 ^ and obtains, 



x = p. 
P = - 



V(x)- I dT7(t-r)p(r)+e(t), 
o 



(4) 



where 7(f) is the memory kernel 
and £(t), the forcing function 

e(*)=X> 



7ft) = ^ — ^ cosw„i, 



Or/fO) ^^(O) ? COS U) v t + 



sin w„t 



(5) 



Having chosen a distribution for the initial state of the bath, given by Eq. © , the fluctuating force £ (t) becomes zero 
centered, and the correlation function of £(£) gives the celebrated fluctuation dissipation relation (FDR) 25 i 26 i 27 



<£(*)£(*)> = fc B T 7 (t-f). 



(6) 



To complete the identification of Eq. |3j) as a generalized Langevin equation, one must establish the conditions on the 
coupling coefficients g v , on the bath frequency w„ and on the number N of the bath oscillators which ensure that 
j(t) is indeed dissipative. Sufficient conditions for 7 (f) to be dissipative are that it is positive definite and decreases 
monotonically with time. Both the conditions are achieved if N — > 00 and if g v /m v Lo1 and u) v are sufficiently smooth 
functions of As N — > 00, one replaces the sum by an integral over uj weighted by a density of states D(lj) to get 



with (g v I 'm u u>~,) — > c(w). For 



7(i) = / duoD(uo)c(uj) cos(wi), 



J D(w)c(w) 



1 + T 2 W 2 ' 



(7) 



(8) 



which can be achieved by a variety of combinations of the density of states D(lu) and the coupling function c(u>), and 
which broadly resembles the behavior of hydrodynamic model in a macroscopic system,— the dissipation kernel 7 (t) 
becomes 



lit) 



7 



exp(-|i|/T c ). 



(9) 



t c in the above expression is the cut-off frequency and is characterized as the correlation time of the bath. In the 
limit t c — > 0, 7(4) — > 27<5(t) and one obtains the traditional Langevin equation in the Markovian domain 



-y (ar)-7P + C(t), 



(10) 



where, (^(i)) = and (£(t)£(t')) — 2 7 fcsT(5(i — t'). If one considers the dynamics of the Brownian particle in aperiodic 
potential V(x) = V(x + L), whose spatial symmetry can be broken by an external load (force) thereby creating a 
biased force field, then the system's dynamics is governed by 



x = p, 



(11) 



where F is the external force. The sum of the periodic potential V(x) and the potential —Fx due to the external 
force F, i.e., U(x) = V(x) — Fx, is a corrugated plane whose average slope (a measurement of the bias) is determined 
by the external force F£2- 

Eq. (| 11 p is the standard Langevin equation of a particle moving in an external potential under an external load force 
and is Markovian in nature. In addition to that, the dissipation term 7 is constant due to the linear system reservoir 
coupling g v and the noise term £(t) is Gaussian, additive in nature reflecting the Markovian kinetics of the Brownian 
particle. In the following subsection, we show how this Markovian kinetics changes to a non-Markovian one due to 
the sudden excitation of the few bath modes and splits the noise term £(t) into two parts. 
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B. The nonstationary system reservoir model 



We consider a Brownian particle of unit mass, described by the coordinate x and the conjugate momentum p, 
moving in a periodic potential of periodicity L, i.e. V(x + L) = V(x). It is acted upon by an external force F, 
which for the present study is assumed to be constant and time independent. The system mode is coupled to a 
set of relaxing modes considered as a semi- infinite dimensional system ( {(^-subsystem) which effectively constitutes 
a nonequilibrium bathi 14 i 19 ' 22 These {g^} modes are in contact with a thermally equilibrated reservoir. Both the 
reservoirs are composed of two sets of harmonic oscillators of unit mass characterized by the frequency sets {ujk} and 
for the nonequilibrium and the equilibrium bath respectively. The system reservoir combination evolves under 
the total Hamiltonian 

H = y + V{x) -Fx+ l - £(P? + fi?Q?) + \ 5>fc + «M - ~ 9(?) E 9k - E a ik9kQ S - (12) 

j k k k 

In Eq. (|12[) . Kj is the coupling constant describing the coupling of the system with the equilibrium bath modes and g(x) 
is the coupling function. The term g(x) J2j Qj indicates the coupling of the nonequilibrium bath to the system and 
the last term describes the coupling between the nonequilibrium bath and the thermal bath with coupling constant 
Ojfc. The equilibrium bath is assumed to be in thermal equilibrium at a temperature T and the initial distribution of 
equilibrium bath variables are assumed to Gaussian. The form of the nonequilibrium bath, that of a set of phonons 
or photons, is chosen for both simplicity and because of its generic relationship to many condensed matter type 
systems j2£ 

Eliminating the equilibrium bath variables {Qj, Pj} in the traditional wa y 25 i 26 ' 27 one may show that the nonequi- 
librium bath modes obey the dynamic equations 

4k = Pk, 

Pk = ~7Pk - ^tlk - g{x) + rj k (t). (13) 

Eq. (|13|) takes into account the average dissipation 7 of the nonequilibrium reservoir modes qt due to their coupling to 
the thermal reservoir which induces fluctuations f]k(t) characterized by the usual FDR (i]k(t)rik(0)} = 2jksT5(t)M^ 
In general, {r)k(t)} being a non-zero constant quantity which, without loss of any generality, may be chosen as zero by 
shifting the origin of our coordinate system as we are dealing with a periodic potential. In passing we mention that 
in deriving Eq. (fT3|) from Eq.(fl"2|). the cross terms for J2j IkjQj have been neglected. 

Proceeding similarly to eliminate the thermal reservoir variables from the equations of motion of the system, we 
obtain 



x = p, 

P = -~feP-V'(x) + F + Z e {t)+g'(x)Y,<lk. (14) 

k 

where 7 e refers to the dissipation coefficient of the system mode due to its direct coupling to the thermal bath 
providing fluctuations £ e (t). The statistical properties of £ e (t) are (£ e (*)) = and (£ e (£)£ e (i )) = 2j e kBT5(t — t'). 
Comparing with Eq. lfTTj) . it is easy to see that the dissipation term j e and the noise term £ e (i) are basically 7 and £ (t), 
respectively, that arise due to the direct linear system reservoir coupling. Now making use of the formal solution of 
Eq. (|13p which takes into account the relaxation of the nonequilibrium modes, and integrating over the nonequilibrium 
bath with a Debye type frequency distribution of the form 19 

= for \u\ > uj c , (15) 

where lo c is the high frequency Debye cut-off, one finally obtains the following Langevin equation for the system mode, 
from Eq. (fT4| as 

x = p, 

p = -T(x)p-V'(x) + F + Ut)+9(^ n {t). (16) 

In the above Eq. (fTrJ)) 

T(x) = le+ln [g'(x)} 2 , (17) 
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is the state dependent dissipation constant comprising of 7„ and 7 e . refers to the fluctuations of the nonequilibrium 
bath modes which effectively cause a damping of the system mode. This damping is also state dependent due to the 
nonlinear coupling function g(x) as is given by j„\g (x)] 2 . In Eq. (fT6|) . the potential V(x) in which the particle moves 
has been modified to 



V{x) = V(x) - ^ ln g\x). 



(18) 



The fluctuations £ n (i) due to the presence of nonequilibrium bath is also assumed to be Gaussian with zero mean 
{£"(*)) — 0- Also, the essential properties of £, n (t) explicitly depend on the nonequilibrium state of the intermediate 
oscillator modes {qj} through u(u,t), the energy density distribution function at time t in terms of the following FDR 
for the nonequilibrium bathi 9 - 



»(w,t) = ~r~ f dT(Ut)Ut- 

47n J-oo 



= -k B T 
2 



-ft/2 



u(u,Q) - -jk B T 



(19) 



[u{uj, 0) — (fcgT/2)] is a measure of the departure of energy density from thermal average at t = 0. The exponential 
term exp(— jt/2) implies that this deviation, due to the initial excitation, decays asymptotically to zero as t — > oo, 
so that one recovers the usual FDR for the thermal bath. 19 ' 22 Eq. (TT9")) thus attributes the nonstationary character of 
the {(?fc}-subsystem. At this point it is pertinent to note that the above derivation is based on the assumption that 
£ n (t) is effectively stationary on the fast correlation time scale of the equilibrium bath modes. This is necessary for 
the systematic separation of the time scales involved in the dynamics. 

Eq. (|16j) is the required Langevin equation for the particle moving in a modified potential V{x) and is acted upon 
by a uniform force F. The motion of the particle is damped by a state dependent friction T(x). Depending on the 
coupling function g(x), both V(x) and r(x) are, in general, nonlinear in nature. As a result, the stochastic differential 
Eg. (|16[) becomes nonlinear. The fluctuating part in Eq. (|16p is comprised of two quantities; t; e (t), an additive noise 
due to thermal bath and £ n (i)> & multiplicative noise due to nonlinear coupling to the {qfc}-subsystcm. The Langevin 
equation (|16p describes a non-Markovian process as well, where the non-Markovian nature is characterized by the 
decaying term in Eq. (JTHJ) , describing the initial nonequilibrium nature of the {gfc}-subsystem created by applying 
sudden excitation at t = . 19 ' 22 



III. STOCHASTIC DYNAMICS IN THE OVERDAMPED REGIME AND THE TIME DEPENDENT 

DISTRIBUTION 



For large dissipation, i.e., in the overdamped limit one usually eliminates the fast variable p adiabatically by 
omitting the inertial term dp/ dt from the dynamical equations of motion to get a simpler description of the system 
in position space. The approach of adiabatically eliminating fast variables is valid on a much slower time scale and 
is a zero order approximation. For constant large dissipation, this adiabatic elimination of the fast variables leads to 
the correct description of the system's dynamics. However, in presence of hydrodynamic interactions, i.e., when the 
dissipation is state dependent, the traditional adiabatic reduction of fast variables does not work properly and gives an 
incorrect description of the system's dynamics. For state dependent dissipation, an alternative approach was proposed 
m RefH Using the method given in ReflU, and using Eq.(TTS]) one may carry out a systematic expansion of the 
relevant variable in powers of l/7 e by neglecting terms smaller than 0(l/~f e ). Then, by Stratonovich interpretation, it 
is possible to obtain the ap pro priate Langevin equation corresponding to a Fokker-Planck equation (FPE) in position 
space. Thus, following Refl24l. the formal FPE for the probability density function (PDF) P(x, t) corresponding to 
the process described by Eq. fTBI) can be obtained as 

dP d (V'{x)-F \ ,.. Lrr df 1 d 1 d1 i 7„ ( t i ^t/2\ d \g(x) d 



InkBT 1 1 + re 




dt dx [ T(x) J ' dx {T(x) dxT(x) J '""V J dx [T(x) dx T(x) 

where r = {[u(ui — > 0, Q)/2ksT] — 1} and is a measure of the deviation from equilibrium at t — 0. Under the steady 
state condition (at t — + oo), dP/dt = and the stationary distribution obeys the following relation, 

fcsT^M + (v'(x) F) P s (x) = 0, (21) 
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which has the solution 



P s (x) = TVexp 



(V'(x) 



F)dx 



(22) 



where N is the normalization constant. In Stratonovich description, the Langevin equation corresponding to the FPE 
given by Eq. (|2Tj) is 



(V'(x)-F) D(t)g'(x)g"(x) 



T(x) 



T 2 (x) 



1 . . . o (x) , , . 
&(*) + T^«(*), 



r(x) 



r(x) 



(23) 



with D = ^f n k B T (1 + r exp(— ^t/2)) being the time dependent diffusion constant due to the relaxation of nonequilib- 
rium bath modes.— Let us consider that the time dependent solution of Eq. ([20|) is given by^ 



P(x,t)=P s (x)eM-<l>(t)), 
where is a function of time only and lim^oo <f>(t) = 0. Ps{ x ) is the steady state solution of Eq. ([2H)l 



(24) 



d \ (V (x) - F) \ , _ d 

Tx \ T(x) PS{X) ] + ^ BT dx \ T(.r) ,i, r,,. 

. , d J g'{x)g"{x) \ 



1 d 1 . \ d J g (x) d g (x) 



(25) 



Substitution of Eq. ([M|) in Eq. (f2"0|) separates the space and time parts and we have the dynamic equation for <f>(t) 



-— exp(jt/2) — const = a (say). 



On integration over time we get, 



2a 

1 



exp(-7t/2), 



where a can be determined from the initial condition. The time dependent solution of Eq. (|20[) thus reads as 

2a 



P(x,t)=P s {x)e*v 
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■ exp(- 7 t/2) 



(26) 



(27) 



To determine a, we now demand that just at the moment the system (and the non-thermal bath) is subjected 
to external excitation at t = 0, the distribution must coincide with the usual Boltzmann distribution where the 
energy term in the Boltzmann factor, in addition to the usual kinetic and potential energy terms, contains the initial 
fluctuation of energy density Au[= u(u>, 0) — (/cbT/2)]. This demands that 



2k B T 1 



(28) 



a is thus determined in terms of relaxing mode parameters and fluctuations of the energy density distribution at 
i = 0. 



IV. STATIONARY AND TRANSIENT CURRENT 



In the over damped limit, the stationary current from Eq. ([25|) can be represented as 



1 



T(x) 



V {x)-F + k B T 



dx 



Ps(x). 



(29) 



Integrating Eq. (|29p we have the expression for stationary probability distribution in terms of stationary current as 



- P -U(x) 



h(x) 



P S (0) Jsle 



h{0) k B T 



h{x')e u ^'Ux 



(30) 
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where h(x) = 1 + {j n /j e )[g (x)] 2 , T(x) = 7 e /i(x) and U{x) — 7 e / dx'h(x )[V (x ) — F]/k B T. We now consider a 
symmetric periodic potential with periodicity L, i.e. V(x) = V(x + L) as well as the periodic derivative of coupling 
function with the same periodicity as that of the potential, i.e., g (x) = g (x + L). As a consequence of this choice, 
U(x) is also a periodic function of x with the period L. If we impose the condition that Ps(x) is bounded for large 
enough x, it follows from the above mentioned conditions of periodicity, that Ps(x + L) = Ps(x) i.e. Ps(x) must be 
periodic with the same period Now applying the periodicity condition of Ps(x), we have from Eq. (j3"P)) 



Ps(0) 

MO) 



~f e /k B T f u{x) 
J s l _ e u(L) J o ( >dx - 



Because of the periodicity, we normalize the steady state PDF in the periodic interval 

f-L 



f P s {x)dx = 1, 
Jo 



to get 



h{x)e 



-U{x) 



Ps(0) Jsle 



h(x')e u{x,) dx 



dx — 1. 



h(0) k B T Jo 

Now eliminating Ps(0)/h(0) from Eq. (j3Tj) and Eq.([33|). one obtains the steady state current 



(31) 



(32) 



(33) 



Js 



h(x)e- u{x) dx / h(x')e u ( x 'Ux 



- e u{L) 



h(x)e- u{x) / h(x')e u(x ' ) dx')dx 



(34) 



From Eq. ((34|) it is clear that in the absence of any external bias F, the steady current vanishes. We thus observe that 
there is no occurrence of current for a periodic potential and for periodic derivative of the coupling function with the 
same periodicity for F = 0. At the macroscopic level this confirms that there is no generation of perpetual motion 
of the second kind, i.e., no violation of second law of thermodynamics. In passing, we note that in the absence of 
{<7/c}-subsystem, i.e., when 7„ = , Eq. ([34"f reduces to the standard forrn^ 



J s = L-f e k B T 



1 _ e LF/k B T 



Mx)/k B T dx 



-V(x)/k B T dx _U_ £ -2LF/k B T 



V( X )/k B T / e V ix )/k B T dx dx 



(35) 



Next, to find the time dependent current J(x,t) we resort to Eq. ([20|) and observe that 



d J V(x') - F 



T(x) 



P } + le k B T 



1 d 1 



r(a:) dxT{x) 



P > + lnk B T 1 + re 



-ft/2 



g'(x) d g'(x) 
T{x) dx T(x) 



+^f n k B T (l 



- 7 t/2\ j 9{x)g"{x) 

r 2 (x) 



re 



P 



(36) 



Now substituting Eq. (|24"|) in Eq. ((36| and making use of Eq.(|25|) we find that J(x, t) can be expressed in a much simpler 
form 



J(x,t) = J s e-W> - V{t)-±-±-±-[g {x)] 2 P s {x), 

I (x) dx I (x) 

where Ps{x) is the steady state PDF and Js is the steady state current given by Eq. 



and 



(37) 



(38) 
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The steady state current Js thus can be obtained from 



(V (x)-F) 1 d 1 . g (x) d g (x) . g {x)g (x) 

J s = ^7-^ Ps(x) - lek B T—— — ——P s {x) - lnk B T—— — ——Ps{x) - "f n k B T — — — — P s {x),{39) 



T(x) 

from which we have 



1 d 1 



■\g(x)rP 8 (x) 



T(x) dx T(x) r(x) dx T(x) 

Js (V'(x)-F) 7e 1 d 1 



r 2 (x) 



r(x) dx T(x 
Using Eq. ([40|) we then obtain from Eq. ([37l 



j n k B T j n k B TT(x) 7„ r(ar) dx T(x) 



Ps(x). 



J(x, t) = J s 



lnk B T 



Jnk B T 



{V ' {x) ~ F Wx) + , e k B T^±^Ps { x) 



T(x) 



T(x) dxT(x)' 



Defining the space dependent part on the RHS of Eq. (|4"Tj) as M(x), we obtain 



J{x,t) = J s 



,-<t>{t) 



l n k B T 



lnk B T 



(40) 



(41) 



(42) 



where 



M(x) 



r(x) 



T(x) dxT(x)' 



(43) 



From Eq. (|4"2")) we observe that the current J(x,t) can be written as a sum of two terms. The first term is space 
independent and only a function of time. The second term is product separable in the form of time and space part. 
As t — > 00, RHS of T>(t) — > and asymptotically J(x,t) reduces to the steady state current Js- Now using the 
continuity equation 

dP(x,t) dJ(x,t) 



dt 

along with P(x, t) = P s (a;)e~* (t) , we get from Eq.® 

dM(x) 



a 



dx 



dx 



Ps(x), 



(44) 



or equivalently 



M{x) 



Ps{x)dx. 



(45) 



As we are dealing with periodic functions, the constant of integration is chosen to be zero. Now integrating Eq. 
for Ps(x) and using the normalization condition, Eq. (j32|) . we have the expression for steady state PDF as 



p s ( x ) = e -(V(x)-F*)/k B T 



1 + M_ jL e -(V(x)-F X )/k B T { J* T ( x ') e (V(x')-F X >)/k B T dx > }dx 



(46) 



/o e-( v W- F *y k BTdx 

Using Eq. p5|) along with Eq. (|4"S")) one obtains from Eq. (|4"2")l the expression for the time dependent current, J(x, t) as 



J(x,t) = J s 



ae 



- 7 t/2 J x dx ' e -(V(x')~Fx')/k E 



^ e -(V(x)-Fx)/k B T dx 

x<j / T(x"){e {v ^^- Fx '"y k - T dx"}dx' 



Js 
k B T Jo 



-{V[x )-Fx )/k B T 



(47) 
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where Js is given by Eq. (l34"l) . Since the potential possesses spatial periodicity, one has J(x,t) = J(x + L.t). Hence 
the net time dependent current is given by 



L jo 



j(*) = T / J{x,t)dx. (48) 



It should be noted that for symmetric potential with F = 0, Js = 0. But, in our development, transient current exists 
and the direction of current depends on the sign of a. What is immediately apparent is that for symmetric potential, 
the sign of Aw[= u(ui, 0) — (l/2fc,gT)] determines the direction of initial current 

j( t) = ae ^ t/2 £ dx $ X dx ' zM-V{x)/kBT] 
L dxexp[-V(x)/k B T) 

It is also clear from Eq. (l47p that the time dependent current reduces to the steady state current, Js in the asymptotic 
limit. The presence of the term exp[— 4>{t)\ in the expression of J(x,t) makes the transient current strongly non- 
exponential in nature. The transient behavior of growth or decay of charge and current in L — R, C — R or L — C — R 
circuit is important in construction of many electrical and electronic devices where there is the mechanism of storage 
of energy. In construction of molecular motor or nano-switch, the transient behavior of the devices may be worth 
studying. In our development, the preparation of intermediate relaxing bath plays a key role to generate the time 
dependent current. Nevertheless, our methodology will also be applicable in the case when any arbitrarily prepared 
bath is approaching towards equilibrium. In passing, we mention that the model considered in the present paper may 
be realized in a guest-host system embedded in a lattice where the immediate neighborhood of the guest comprises 
intermediate oscillatory modes, while the lattice acts as a thermal bath. 



V. CONCLUSION 



We have hereby proposed a simple microscopic system nonequilibrium bath model to simulate nonstationary 
Langevin dynamics. The nonequilibrium bath is effectively realized in terms of a semi-infinite dimensional reser- 
voir which is subsequently kept in contact with a thermal reservoir which allows the non-thermal bath to relax with 
characteristic time. The frequency spectrum of the relaxing bath is assumed to be of Debye type. By an appro- 
priate separation of time scale, we then construct the Langevin equation for a particle in which the dissipation is 
state dependent and the stochastic forces appearing are both additive and multiplicative. The underlying stochastic 
dynamics is found to be nonstationary and non-Markovian. Based on the strategy of Sancho et al.,-- we then show 
that this Langevin equation can be recast into the form of generalized nonstationary Smoluchowski equation which 
reduces to its standard form asymptotically. We then solve the expression for time dependent PDF. As an immediate 
application of our recent development, we consider the dynamics of a Langevin particle in a ratchet potential and 
obtain the analytic expression for both stationary and nonstationary transient average velocity, which is followed by an 
immediate observation that in a periodic potential the direction of nonstationary current depends on the preparation 
of the nonequilibrium bath. 
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